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Abstract. Spherical meteoroids orbiting the Sun experience orbital decay due to the Poynting- 
Robertson effect and can become trapped in commensurability resonances with a planet. Detail 
conditions for trapping in the planar circular restricted three-body problem with action of solar 
electromagnetic radiation on spherical grain are presented. The results are given in terms of the 
radiation presure factor f3. 

In general, the greater value of /?, the shorter capture time in a resonance. Captures of 
particles may be even several times longer than it is presented in literature. Analytical formula 
for maximum capture times is given. 

Analytical theory for minimal capture eccentricities is presented. The obtained analytical 
results differ from those presented in the literature. However, the known analytical results are 
not very consistent with our detail numerical computational results for j3— values 0.01 and 0.05 
and for the planet Earth. Minimum eccentricities for captures into resonances are also smaller, 
even in order of magnitude, than the published values. 

The presented results enable comparison with calculations for nonspherical particles. 
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1. Introduction 

Physics of resonances with Solar System planets is discussed mainly during the last three decades. The 
orbital evolution of meteoroids near resonances with planets was also investigated. Besides gravita- 
tional forces the effect of solar electromagnetic radiation in the form of the Poynting-Robertson (P-R) 
effect is usually taken into account (Robertson 1937, Klacka 2004, 2008a, 2008b), see e. g., Jackson 
and Zook (1992), Weidenschilling and Jackson (1993), Marzari and Vanzani (1994), Sidlichovsky and 
Nesvorny (1994), Liou and Zook (1995), Liou et al. (1995). 

The aim of our contribution is detail analytical reconsideration of the influence of the P-R effect on a 
meteoroid which can become trapped in commensurability resonances with a planet. As a motivation 
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we can introduce that Weidenschilling and Jackson (1993) yields for theoretical minimal value of 
eccentricity at the beginning of capture the value 0.3506 for dust particle with /3 = 0.01 in resonance 2:1 
with the Earth, while our simulations show that real values may be even less than 0.1. Detail conditions 
for trapping in the planar circular restricted three-body problem with action of solar electromagnetic 
radiation on a spherical grain (P-R effect holds) are presented. We concentrate on real values of orbital 
elements - width of semimajor axis in resonance, minimal values of eccentricities at the beginning of 
capture, and, capture time. We use for calculating osculating orbital elements simultaneously solar 
gravity and dominant part of solar radiation pressure-term (central acceleration have form — GM Q (1 — 
[3)ep/r 2 ). Results presented in this contribution may serve for comparison of the results obtained for 
spherical and nonspherical dust particles (meteoroids), since P-R effect holds only if a special condition 
is fulfilled (see Eq. (9) in Klacka 2008a; compare Eq. (73) in Klacka 2008b) and real nonspherical 
particles do not fulfill the condition. 



2. Model - Equation of motion 

All our numerical simulations are based on equation of motion written in the form (see Eq. (72) in 
Klacka 2008b for more details) 

dv GM Q ( r — r P r P 

ep — Gm 



dt r 2 \ \r — rp\ 3 \rp\ 3 

GM P 



A' |m 2 l 

e R = r/\r\ ; /3 = 7.682 x 10~ 4 Q' \— ± , (1) 



where r is position vector of the particle (with respect to the Sun) and rp is position vector of 
the planet which moves in circular orbit around the Sun; mp is mass of the planet, M© is mass of 
the Sun, m is mass of the particle/meteoroid, A' = ir a' 2 , a' is radius of the spherical particle, Q' pr is 
climensionlcss efficiency factor (see, e. g., van de Hulst 1981). j3 = 5.761 x 10" 5 Q' pr / (e>[g/cni ] s[cm]), 
for homogeneous spherical particle in the Solar System, g is mass density and s is radius of the sphere. 
Planar motion will be considered. 

We use osculating orbital elements which considers solar gravity and dominant part of solar radi- 
ation pressure-term, simultaneously: central acceleration — GM Q (1 — (3) ep/r 2 is used in calculations 
of orbital elements and the corresponding quantities will be denoted with a subscript (3. 
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3. Perturbation equations of celestial mechanics 
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Rewriting Eq. (1) into the form 



dv 
~dt 



GM & (1 - (3) 



- P 



GM G 



(- 



e R - Gmp 



r P 



rp\ 



rp 

\r P \* 



ep 



(2) 



we can immediately write for perturbation acceleration to Keplerian motion 
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Perturbation equations of celestial mechanics yield for osculating orbital elements (ap - semimajor 
axis; ep - eccentricity; ip - inclination (of the orbital plane to the reference frame); flp - longitude of 
the ascending node; cop - argument of pericenter; Qp is the position angle of the particle on the orbit, 
when measured from the ascending node in the direction of the particle's motion, <dp = ujp + fp): 
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where \i = GMq and r = pp/(l + ep cos /^) ; Fp p, Fp t and Fp jy are radial, transversal and normal 
components of perturbation acceleration. 
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3.1. P-R effect and perturbation equations of celestial mechanics 



This subsection follows considerations presented in Klacka (2004 - Sec. 6.1), or, in Klacka (1992). 
On the basis of Eq. (3), we can immediately write for components of perturbation acceleration to 
Kcplerian motion: 

{F p R ) pR = -2(3^ V -^f, (F 0T ) pR = -p£ V -ef, (Fp n) pr = , (5) 
where fi = GM Q , r — pp/(l + ep cos fp) and two-body problem yields 
vp r = \J li (1 - P)/pp ep sin fp , 

vpT = (1 - P)/pp (l + epcosfp) . (6) 

Inserting Eqs. (5) - (6) into Eq. (4), one easily obtains 

[i 2ap 1 + e| + 2ep cos fp + e| sin 2 fp 
r z c 1 — ep 

~ P \ (2e ( 3 + e (3 sin 2 / /3 + 2cos/ / 3) , 
, 
, 

- (3 4 - — " (2 - ep cos fp) sin fp , 

r z c ep 

y / ^{l-f3)pp ^ 
3.2. Gravity of a planet and perturbation equations of celestial mechanics 

As for gravitational acceleration generated by a planet (Fp) G presented in Eq. (3), we can write 
(compare Eqs. (44) in Klacka 2004 for the case ip = Qp = ujp = and Eqs. (57) in Klacka 2004) 

r P = r P e PR , r = r e R , 

e PR = (cos9p,sin9p,0) , 

e R = (cos 9,3, sin 9/3,0) , 

ep = (— sin 9/3, cos 9/3,0) , 

ep R = cos (9p - 9/3) e R + sin (9p - Qp) e T , 

9p = n P (t - t ) + 9pq , n P = \J G (M Q + m P ) r p 3/2 . (8) 
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Inserting Eqs. (8) into (Fp) G presented in Eq. (3), we obtain 

(Fp) G = Gm P (Xj - X„) , 

x _ [rp cos(9 P - Qp) - r] e R + r P sin (Op - Op) e T 
1 [rp + r 2 - 2 r r P cos (6 P - Op)} 3/2 

X H = \ [cos(9 P - 0^) ep + sin(e P - @p) e T ] . (9) 
r p 

In order to obtain (dap/dt)o, (dep/dt)a, ■■■ (dQp/dt)a, it is sufficient to put Eqs. (9) into Eqs. (4). 
4. Exterior mean motion resonances and dap/dt 

We will consider only terms proportional to e° and ep, i. e., we will neglect ep and higher orders in 
ep. First of Eqs. (4) can be written in the form 

= 2 0/3 ^{l- J) {Fp R 6/3 Sin//3 + F 0t{^ + e cos/^)} . (10) 
On the basis of Eqs. (9) we write, then (fip = Gmp): 

(r P cos(9 P - 6/3) - ap cos(9 P - Qp)\ 
{Ff3R) ° = MP \ Yjo) r% J 

( F P t) g = Hp 

x (ep) = [r% + ap (1 - 2ep cos fp) - 2apr P (1 - ep cos fp) cos (6 P - Gp)] , (11) 
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for radial and transversal components of gravitational perturbation acceleration. 

A particle is in exterior mean motion resonance with a planet when the ratio of their mean motions 
is approximately the ratio of two small natural numbers: np/np = j j (j + q) corresponds to q— order 
exterior resonance. For the first-order resonances we have np/np = j/(j + 1), i. e., any first-order 
resonance may be defined by the ratio T/T P = (j + 

As for commcnsurability resonances with a planet moving in the circular orbit, the third Kepler's 
law [apnp = u(l - apnp = u(l + mp/M Q )] yields 

^=(i-/?) i/s (ix /3 fi + ^v 1/3 , d2) 



a P K ' \T P J \ M Q/ 

for the semimajor axis and period of revolution around the Sun of the particle and the planet (subscript 
P). Defining any resonance by the ratio T/Tp, Eq. (12) yields immediately the ratio ap/ap. 
On the basis of Eqs. (11) we can write for the important resonant terms 

RES = {{Fp R ) G ep sin fp + (Fp T ) G (1 + ep cos fp)} , 

RES = hp | -^j^ep sin (9 P - up) - ^j^ep sin fp + ^ sin (9 P - 6^ )| . (13) 
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Using the relation 

VI \ ll Q t 1 - £ COS (9p - 6/3) 

XM = X{0) i 1 - 3< * COfl ><' l-2 £ cos(ep-9 /3 )+ £2 
e = rp/ap , 

We can Eq. (13) rewrite to the following form 
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- sin (9^ — 26p + w^) + - sin (3Qp — 26p — up) — 2 sin (Qp — up) 



+ e 3 sin (9p — up) 



(15) 



In this expression we neglect terms proportional to sin(9/3 — 6p), because they have not large influence 
when we take into account terms important in first-order mean motion resonance. We use Laplace 
coefficients bj defined by the relation 

(l-2ecos(G fj -G P )+e 2 )- 5 / 2 = - £ b 3 cos j(Qp - Q P ). (16) 



From Eq. (15) using Eq. (16) we get for the first-order exterior resonance 
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where we have taken into account only terms with sine argument (7 + 1)9/3 — j9p — up — T. We 



define Cj as 



a = - 



.7 + 1 



~2 J U J+1 ~ 4 j_1 
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Using Eqs. (10), (13), (17) and (18) we obtain 
dt 



dap\ MP ~ . „ 

- -J =~ U + !) a /3 e /3 n /3 ^(1-/3) ^ smr 



(18) 



(19) 



For the first-order exterior mean motion resonance Eq. (12) reduces to 
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Table 1. Values of Cj for particle with j3 = 0.01 and (3 = 0.05 in first-order exterior mean motion 
resonance with the Earth in a circular orbit. 



j 


C7 i (/3 = 0.01) 


<7,(/J = 0.05) 


1 


1.31819 


1.35328 


2 


1.68270 


1.77804 


3 


2.05585 


2.23845 


4 


2.43605 


2.73560 


5 


2.82257 


3.27183 


6 


3.21514 


3.85027 


7 


3.61366 


4.47467 


8 


4.01813 


5.14933 


9 


4.42856 


5.87920 


10 


4.84502 


6.66992 


11 


5.26757 


7.52797 


12 


5.69630 


8.46078 


13 


6.13128 


9.47693 


14 


6.57262 


10.58637 



The values of Cj are given in Tab. 1 for particles with (i € {0.01,0.05} in the first-order exterior 
resonance j € {1, 2, 3, 14} with planet Earth in circular orbit. 

Neglecting higher orders of eccentricity e/j, the first of Eqs. (7) yields 



(22) 

PR 



^) =-2 0-0- (1 + 4e /3 cos/ /3 ) . (21) 
at ) PR ape 

Eqs. (3) and (4) give 

dap _ f dap \ ( dap_ 

dt \ dt ) q V dt 

Exterior resonances can produce long-lived trapping: dap/dt = 0. On the basis of Eqs. (19), (21) and 
(22) we can write 

- (j + 1) a fj en n & ^ P C, sinT = 2 p (1 + 4e p cos f p ) . (23) 

Relation np — np + 1) and the last one of Eqs. (8) yield np = y/JI (1 + mp/M Q ) 1 / 2 a p 3 ^ 2 
+ 1). Using this result and also Eq. (20) in Eq. (23), 
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Eq. (24) immediately yields criterion for resonant trapping: 

Sa 4. 1 ^ c a+mp/M Q r i/6 (j + if 3 r 1 

e@ mm = < 4 + -r-p: Cj > . (25) 

5. Comparison of analytical and numerical results 

We have numerically solved Eq. (1) for various values of j3 and for various initial conditions. The most 
important results are presented in Tables 2 and 3 for the case of the first-order exterior resonances 
with the Earth. Numerical results are compared with analytical results obtained from the theory 
presented in Sees. 2-4 and, also, with the analytical results presented by Weidenschilling and Jackson 
(1993). Our theoretical values of semimajor axis for the exterior resonance (j + are obtained 
from Eq. (20). The values are given in columns labeled as ap. The values ap m „ and ap max are 
minimal and maximal semimajor axes during the capture in a given resonance. The width of the 
resonances in semimajor axis can be calculated from the relation ap max — ap m in- Minimal values of 
eccentricities for the capture into resonances, found by numerical solution of Eq. (1), are given in the 
column ep m i n num. - Minimal capture eccentricities ep m i n th. were computed from Eq. (25). Minimal 
capture eccentricities ep m j„ t h. wji presented in Tab. 2, are taken from Table II in Weidenschilling 
and Jackson (1993). Maximal capture times for different first-order resonances are also given. The 
times and the minimum and maximum semimajor axis values are obtained for a given numerical 
integration of a mctcoroid. The extremal values of semimajor axis occur when the maximal capture 
time exists, for a given resonance. The minimum eccentricity values were found independently There 
is no relation between the minimum eccentricity and the set {maximal capture times, minimum and 
maximum semimajor axis values}, as for the presentations in Tables 2 and 3. 

General trend seen from Tables 2 and 3 is that the larger value of j - type of the first-order exterior 
resonances with the Earth is (j + l)/j, the smaller width of resonances in semimajor axes and the 
shorter capture times. Moreover, the greater value of /?, i) the smaller resonant width in semimajor 
axes, and, ii) the shorter capture time. 

6. Extrapolation of capture times 

We found, from numerical integrations of Eq. (1), that maximal capture time for a particle in a given 
mean motion resonance is proportional to the value of the particle's (3 and semimajor axis ap according 
to 

al 

r~f. (26) 
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Table 2. Characteristics for the first-order + resonances with Earth for particle with (i = 0.01. 



3 


ag 


Q>(3 min 


Qf3 max 


&/3 min 


min 


C/3 min 


maximal 










num. 


th. 


th. WJ 


capture time 




[AU] 


[AU] 


[AU] 








[ 10 3 years ] 


1 


1.5821 


1.5761 


1.5881 


0.102 


0.1107 


0.3506 


1492 


2 


1.306 


1.3023 


1.3097 


0.036 


0.0784 


0.0387 


736 


3 


1.2074 


1.2046 


1.2101 


0.011 


0.0564 


0.0227 


742 


4 


1.1565 


1.1544 


1.1586 


0.008 


0.0418 


0.0148 


429 


5 


1.1255 


1.1238 


1.1272 





0.032 


0.0104 


268 


6 


1.1045 


1.1031 


1.1059 


0.003 


0.025 


0.0077 


252 


7 


1.0894 


1.0881 


1.0907 


0.001 


0.0201 


0.0059 


227 


8 


1.0781 


1.077 


1.079 





0.0164 


0.0046 


362 


9 


1.0692 


1.0681 


1.0701 


0.001 


0.0136 


0.0038 


167 


10 


1.062 


1.061 


1.0626 





0.0114 


0.0031 


469 


11 


1.0562 


1.0551 


1.0566 


0.002 


0.0097 


0.0026 


75 


12 


1.0513 


1.0503 


1.0516 


0.001 


0.0084 


0.0022 


36 


13 


1.0471 


1.0459 


1.0479 


0.001 


0.0073 


0.0019 


21 


14 


1.0436 


1.0423 


1.0444 


0.001 


0.0064 


0.0016 


13 



This relation is also in accordance with the expected dependence of secular time derivation of eccen- 
tricity (Liou and Zook 1997; Klacka et al. 2008 - r\ = 0) on ap and j3. Using Eq. (12) (with assumption 
mp -C Mq) we obtain from Eq. (26) 

r(api,/3i) _ (a P i\ 2 /3 2 (l-/3i) 2/3 (27) 
r(a P2 ,/3 2 ) \a P2 ) fa (1 - /3 2 ) 2 / 3 ' 1 ' 

Eq. (27) presents the analytical formula for obtaining maximum capture time r(ap2,P2) for the 
particle of /3 2 and planetary orbital radius ap 2 for a given type of resonance, if maximum capture 
time 7"(api,/?i) is known for for the particle of (3i and planetary orbital radius api for the same type 
of resonance. 

Using Eq. (27) we have calculated the values of capture times in Tab. 4 from values of capture 
times in Tables 2 and 3. We use the values in Tab. 2 in calculations of the expected values in Tab. 3, 
and vice versa. 
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Table 3. Characteristics for the first-order + resonances with Earth for particle with = 0.05. 
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ap 


&(3 min 


&/3 max 


6/3 min 


6/3 min 


maximal 










num. 


th. 


capture time 




[AU] 


[AU] 


[AU] 






[ 10 3 years ] 


1 


1.5605 


1.5555 


1.5656 


0.196 


0.1981 


216 


2 


1.2882 


1.2851 


1.2911 


0.096 


0.1701 


168 


3 


1.1909 


1.1887 


1.1931 


0.061 


0.1422 


187 


4 


1.1407 


1.1389 


1.1426 


0.045 


0.1172 


80 


5 


1.1101 


1.1088 


1.1114 


0.028 


0.096 


56 


6 


1.0894 


1.0879 


1.091 


0.018 


0.0786 


102 


7 


1.0746 


1.0735 


1.0756 


0.016 


0.0645 


62 


8 


1.0633 


1.0626 


1.0643 


0.01 


0.0531 


23 


9 


1.0546 


1.0538 


1.0552 


0.007 


0.044 


20 


10 


1.0475 


1.0466 


1.0481 


0.005 


0.0366 


13 


11 


1.0418 


1.0408 


1.0424 


0.005 


0.0307 


6 


12 


1.0369 


1.036 


1.0376 


0.004 


0.0258 


3 


13 


1.0328 


1.0317 


1.0334 


0.005 


0.0218 


5 


14 


1.0293 


1.028 


1.0301 





0.0185 


1 



7. Conclusion 

The planar circular restricted three-body problem with action of solar electromagnetic radiation on a 
spherical grain is considered. This Poynting- Robertson effect is studied numerically and analytically. 
We have found minimal values of capture eccentricities into the first-order exterior mean motion 
resonances. Unfortunately, our theory for the minimal values of the capture eccentricities cannot 
explain the found numerically values. The numerical values ep m i n are also smaller than the values 
presented in Weidenschilling and Jackson (1993). Thus, any analytical method explaining minimal 
values of capture eccentricities does not exist, for the first-order exterior mean motion resonances. 

Our numerical results offer several resonant characteristics for 0— values 0.01 and 0.05 for the first- 
order exterior resonances with the Earth; capture times for = 0.1 are in order of magnitude shorter 
than for = 0.05. We have succeeded in finding greater values of capture times than it is presented 
in Sidlichovsky and Nesvorny (1994). Eq. (27) presents an analytical formula for obtaining maximum 
capture time for the particle of 02 and planetary orbital radius ap2 for a given type of resonance, if 
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Table 4. Extrapolation of capture times presented in Tables 2 and 3. Eq. (27) is used. 



3 


extrapolated 


extrapolated 




values in Tab. 2 


values in Tab. 3 




[ 10 3 years ] 


[ 10 3 years ] 


1 


1110 


290 


2 


863 


143 


3 


961 


144 


4 


411 


83 


5 


288 


52 


6 


524 


49 


7 


319 


44 


8 


118 


70 


9 


103 


32 


10 


67 


91 


11 


31 


15 


12 


15 


7 


13 


26 


4 


14 


5 


3 



maximum capture time is known for for the particle of 0i and planetary orbital radius api for the 
same type of resonance. 

The presented values can be useful in comparison of resonant evolution for spherical and non- 
spherical meteoroids. Our results, based on the P-R effect, can be used for spherical particles, only. 
Real nonspherical particles may exhibit nonzero values for other two components of the radiation 
pressure efficiency factor vector Q' 2 and Q' 3 (Klacka 2004), and, moreover, values of /?— parameter are 
time dependent. Thus, the change of semimajor axis ap in orbital resonance caused by the change 
of parameter (3 may be relatively large. Since the resonant width is a decreasing function of (3, one 
should await that resonant capture times for nonspherical particles are smaller than the capture times 
for spherical grains. 
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